Demonstration of Deutsch's Algorithm on a Stable Linear-Optical Quantum Computer 
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system. By employing photon polarization and spatial modes, we implement all balanced and 
constant functions for a quantum computer. The experimental system is very stable, and the 
experimental data are in excellent accordance with the theoretical results. 

PACS numbers: 03.67.Lx, 03.67.Mn, 42.50.Dv 



a zhang.pei@stu.xjtu.edu.cn 



2 



Quantum computation may solve some complex computational problems and hit the security of the classical cryp- 
tography. It has attracted much interest to investigate quantum algorithms and to realize quantum hardware, which 
are very important to quantum information processing and quantum computation. The first quantum algorithm was 
proposed by Deutsch in 1985 [l[, then extended by Deutsch and Jozsa in 1992 0). The realizations of quantum 
Deutsch's algorithm on quantum computers have been examined in many physical systems, including ion traps [3j, 
nuclear spins in magnetic resonance jj|, super conducting resonators semiconductor quantum dots @, neutral 
atoms Q, and linear optics [8l-fl(]l|. In linear optical systems, it is easy to deal with entanglement and decoherence, 
and the incorporation of detection and post selection make it possible to achieve all-optical quantum computers 
[llj , so linear optical systems are a good candidate for implementing quantum algorithms [IH, [l3| . The system of 
single-photon few-qubit has been used to build the deterministic quantum information processor ( QIP), a nd few-qubit 
QIPs have drawn much attention for application in quantum optics and quantum computation [14l420| . Oliveira et 
al. experimentally tested Deutsch's algorithm using a single-photon two-qubit (SPTQ) system in 2005 [9]. However, 
the light source in their experiment was a bright coherent light instead of a single photon, and they realized the four 
relevant operations using a phase-sensitive Mach-Zehnder interferometer, which needs additional phase stabilization. 
In this Brief Report, we experimentally demonstrate Deutsch's algorithm using a more robust setup [2(| at the single- 
photon level. By employing photon polarization and spatial modes as a SPTQ system, we implement all balanced 
and constant functions for quantum computer. The experimental system is very stable, and the experimental data 
are excellent in accordance with the theoretical results. Furthermore, we also introduce a phase variation in the input 
spatial qubit, which helps us easily to study the differences of all input states for the algorithm. 
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FIG. 1. Quantum circuit of Deutsch's algorithm. H is the Hadamard gate. 



Deutsch's algorithm combines quantum parallelism with a property of quantum interference. Suppose we are given 
a boolean function /(x), where x is either or 1. What we want to know is whether f{x) is a constant function or 
a balanced function, where constant function means f(x) — or f(x) — 1 [or /(0) = /(l)] and balanced function 
means f(x) — x or f(x) = inv(x) [or /(0) ^ /(I)] and 'im;' is the inversion operation. The classic computer has to 
run the f(x) twice to distinguish a balanced function from a constant function, while a quantum computer does the 
job in just one go. Fig. 1 is the quantum circuit implementing Deutsch's algorithm [2l[. Uf is the quantum operation 
which takes inputs \x, y) to \x, y f(x)). A brief explanation is given subsequently. The initial state is l^o) = |0) |1). 
After the Hadamard transformation (H), we get = (|0) + |1))(|0) - |l))/2. Applying Uf to we obtain 

to be one of two possible states, depending on f(x): 

f±(|0> + |l»(|0)-|l»/2, /(0) = /(l), 
1 2> \±(|0>-|l»(|0)-|l»/2, /(0) ^/(l). 1 ' 

The final Hadamard gate is applied on the first qubit, 

f±|0)(|0)-|l»/V2, /(0) = /(l), 

1 3; \±|i)(|o)-|i»/V2, /(o)^/(i), 1 ] 

so we can determine f(x) to be balanced or constant by only measuring the first qubit once. 

From the preceding description, to physically test the algorithm, we need a device which can implement the co- 
operations for the four possible functions. All the possible f(x) functions and Uf operations are summarized in Table 
I. In the first case oiUf = /, it means that the second qubit never changes, whether the first qubit is or 1, so this 
can be recognized as an identity operation to the two qubits. The second case shows that Uf is a NOT gate. The 
second qubit always flips, no matter what the first qubit is. In the third case, Uf is a controlled-NOT (CNOT) gate. 
The second qubit flips when the first qubit is 1. In the last case, Uf is a zero-controlled-NOT (Z-CNOT) gate, where 
the second qubit flips when the first qubit is 0. For these four different Uf operations, identity operation and NOT 
operation are very simple to be realized, and the Z-CNOT gate can be obtained from a CNOT gate with some small 
changes. So the CNOT gate is the fundamental and essential part to execute Deutsch's algorithm. In this context, 
we start with a CNOT gate realized by employing polarization and spatial positions of photons [20], construct the 
four different gates and U / operations, and carry out Deutsch's algorithm. 
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TABLE I. Four different cases of Deutsch's algorithm 



Class 

Constant 
Constant 
Balanced 
Balanced 



Function 

fix) = 
fix) = 1 
fix) = x 
f(x) — inv(x) 



Operation 

\x,y) -» \x,y 0} 
\x,y) -> \x,y®l) 
\x,y) -> \x,y®x) 
\x,y) -> \x,y®(x®l)) 



I (Identity) 
NOT 
CNOT 
Z-CNOT 



The CNOT gate is shown in Fig. 2. The dove prism (DP) is inclined at a 45° angle relative to the horizontal 
plane [shown in Fig. 2(a)], so the images which pass through it from left to right will be rotated by 90°. Suppose 
the polarized beam splitter (PBS) here transmits horizontal-polarized (H) photons and reflects vertical-polarized (V) 
ones. So the H photons travel counterclockwise, while the V photons travel clockwise. With a DP inclined at 45°, the 
spatial mode of H (V) photons is oriented 90° (—90°). Specifically, the left-right (l-r) section of the input photons 
is rotated into the down-up (d-u) section of the output beam for H photons but into the u-d section for V photons 
[shown in Fig. 2(b)]. If we define photon polarization as the control qubit (V — > and H — > 1) and spatial mode as 
the target qubit (Z&u — > and rk,d — > 1), the CNOT operation can be described as follow: 

\V)\l)^\V)\u), \ V )\r)^\V)\d), \H)\l)^\H)\d), \H) \r) -> \H) \u) . (3) 

For the Z-CNOT gate, we should realize the following transition: 

\V)\l)^\V)\d), \V)\r)^\V)\u), \H)\l)^\H)\u), \H) \r) -> \H) \d) . (4) 

Similar to the implementation of the CNOT gate, if we set DP at —45°, the spatial mode of H (V) polarized photons 
will be oriented -90° (90°), and it will be a Z-CNOT gate. The CNOT (Z-CNOT) gate is a polarization Sagnac 
interferometer in our setup, and the two counter-propagating photons always undergo the same amount of phase 
disturbance. So this optical CNOT (Z-CNOT) gate has an inherent stability which requires no active stabilization. 




FIG. 2. (Color online) (a) Optical implementation of CNOT gate by employing polarization and spatial positions of photons. 
The dove prism (DP) is inclined at 45°. Red lines (heavy gray) show the left (I) and right (r) spatial modes, and green lines 
(light gray) show the up (u) and down (d) spatial modes, (b) Spatial positions of input and output beams on the splitting 
plane AB of the polarized beam splitter (PBS). 

We experimentally realize Deutsch's algorithm using the CNOT gate mentioned earlier. The experimental setup 
is shown in Fig. 3. Our source is a He-Ne laser (MELLES GRIOT, 05-LHP-I7I) with deep attenuation to about 
150,000 photon counts per second, which means that the mean distance between two photons is about 2000 m (much 
bigger than our experimental setup length 0.5 m), and the two-photon probability is 2.5 x 10~ 4 . All the PBS are 
quasi symmetric and transmit H photons while reflecting V photons. A polarizer and half wave-plate (HWPi) are 
used to prepare photon polarization states. Here we prepare the initial polarization of photons as V. A 50% beam 
splitter (BS) and a mirror (M) are used to prepare the photon spatial-mode states. The piezo-transmitter (PZT) on 
the first mirror is used to control the relative phase tp between two spatial modes. HWP2 and HWP3 at 22.5° arc 
used as the polarization Hadamard gates. The state after HWP2 can be written as: 

|^i> = (|V0 + |fl))(|0+^|r))/2; (5) 

in particular when (p — tt, \tpx) is equal to l^i), which is mentioned erlier. So this single-photon two-qubit state can 
be used as the input state of Deutsch's algorithm as we described in Fig. 1. Then this state will be evolved by the 
Uf operation. The detection part consists of a Hadamard gate (HWP3), PBS2, and two single-photon detectors (Di 
and D2), which detect the photon's polarization state (the first qubit of the output state). The key point to carry out 
Deutsch's algorithm is how to realize the four different cases of Uf operation. We will discuss these four U / operations 
subsequently. 
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FIG. 3. (Color online) Experimental setup of Deutsch's algorithm. The source is a He-Ne laser (MELLES GRIOT, 05-LHP-171) 
of wavelength 632.8 nm. We attenuate the coherent light into a single-photon level by using some neutral attenuators denoted 
by ATT. P denotes the polarizer for initial state preparation. Three half wave-plates (HWP) and two polarized beam splitters 
(PBS) are used in the setup. A PZT actuator is used to modulate the phase <p between I and r paths. The dove plate (DP) is 
set at 45°. Detectors (Di and D2) are single-photon counting modules (SPCM-AQRH-14-FC). All the mirrors are marked as 
M. 



In the constant-function case, U / can be an identity or NOT operation. For an identity operation, we can simply 
remove PBSi in our setup and set DP at —45°. Therefore, photons in I or r will always undergo a counter clockwise 
route and be output in u or d, respectively, without the effect of polarization. This means that the target qubit 
(spatial mode of photons) will not change with control qubit (polarization of photons). We can deduce the process 
as follow: 

-^\^) = (\V) + \H))(\u)+e^\d))/2 

H ^ 3 m = \V)(\u)+e^\d))/V2, (6) 

whereas for the NOT operation, we can remove PBSi in our setup and set DP at 45°. Then |Z) is converted to \u) 
and \r) is converted to \d). Applying a Hadamard gate (HWP3), we can obtain 

|V> 3 > = \V) (\d) + e* \u))/V2. (7) 

For the preceding two cases, we can only detect the polarization qubits; the results are same without any changes 
when we adjust the relating phase tp. So in our setup, when the boolean function f(x) is a constant function, the 
detector D2 will be clicked, and no photons will arrive at Di. Fig. 4(a) shows our experimental results of Uf = I, 
and Fig. 4(b) shows the results of Uf = NOT. Because there is no interference in these processes, the counts of Di 
and D2 do not change while modulating the voltage of PZT. 

In the balanced-function case, f(x) = x or f(x) = inv(x). We need to place the PBSi into the optical route, so 
the V photons and H photons will travel through the DP in different directions. As we have discussed, when we set 
the DP at 45° (-45°), this will be the CNOT (Z-CNOT) gate for the input state \ipi) shown in Eq. (5). Using the 
corresponding relations of Eq. (3), the theoretical analysis oiUf — CNOT is shown below. 

°™ T life) = (|*0 |«) + e* \V) |d) + \H) \d) + \H) \u))/2 

H ^ a 1^3) = [(1 + e*) \V) (\u) + |d)) + (1 - e*) \H) (|«) - |d))]/V2. (8) 

For the Z-CONT operation, we set the DP at —45°. The output state is 

|^3) = [(1 + e'>) \V) (\u) + |d» (1 - e*) \H) (|u) - \d))]/2V2. (9) 

For these two operations, we still detect the polarization qubits. Then we can get two curves which show the photon 
counts of two detectors changing with the relative phase between two spatial modes. Fig 4(c) corresponds to the 
CNOT operation, and Fig. 4(d) corresponds to the Z-CNOT operation. From Eq. (8) and Eq. (9), we know that the 
theoretical results are sinusoidal functions, and our experimental data fit them well. 

Our experimental results are shown in Fig. 4. In our experiment, we make the relative phase tp adjustable by using 
a PZT controller, so the output state \tjj3) contains the phase parameter ip. When using a PBS for the projective 
detection, the detectors of Di and D2 detect photons of different polarization: H on Di and ^ on D2. From the Eq. 
(8) and Eq. (9), we can see that the photon counts of Di and D2 will sinusoidally vary with the ip being continuously 
changed. We set the phase range for two periods (the voltage of PZT is adjusted from to 34 V) and plot the 
counts-voltage curves. From the description of Deutsch's algorithm, the input state is a certain state with certain a 
phase [Eq. 1]. However, we can get this state simply by setting the phase ip = (27V + \)ir (adjust the PZT in proper 
voltages), where N is an integer. Then Eq. (8) and Eq. (9) are changed into \ip3)^ = ± \H) (\u) — \d))/y/2, where 
'+' is for the CNOT operation and '-' is for the Z-CNOT operation. And if we also set <p = (2N + 1)tt in Eq. (6) 
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FIG. 4. (Color online) Experimental data of Deutsch's algorithm. Black squares show the photon counts of Di, and red 
triangles show the photon counts of D2. Fitting lines are also shown. By modulating the voltage of PZT from to 34 V 
and every 1 V as a step, we record the photon counts of Di and D2 simultaneously, (a) Identity operation of Uj for constant 
function f(x) = 0. (b) NOT operation of Uf for constant function f(x) = 1. (c) CNOT operation of Uj for balanced function 
f(x) = x. (d) Z-CNOT operation of Uf for balanced function f(x) = inv(x). Green dashed vertical lines are used to mark the 
proper points (phases) of the initial states, which can be used to perfectly discriminate the two kinds of functions. 



and Eq. (7), we get — ^ 1^0 (\ u ) ~ l^))/v2) where '+' is for the NOT operation and ' — ' is for the / operation. 

These results are the same as those for ^3), described in Eq. (2). These proper points for Deutsch's algorithm are 
marked by green lines in Fig 4. From these points, we can claim that it is a constant function when Di clicks and a 
balanced function when D2 clicks. Our data also show that we can only probabilistically discriminate the function if 
tp 7^ {2N + l)7r; In particular, when ip — 2Nn, we cannot discriminate the two kinds of f(x) at all. 

Benefiting from the Sagnac interferometer, our experimental setup is very stable without any other additional 
feedback control. This long time stability makes it possible to change the voltage 1 V as a step from to 32 V. We 
can define 77 = I Cnl+Col I as a con trast ratio to describe the precision of our results, where Cui and C02 denote to 
the photon counts of Di and D2, respectively. Theoretically, the contrast ratio r\ is equal to 1. In our experiment, 
for the constant functions, r) c = 99.96 ± 0.03% in Fig. 4(a) and r) c = 99.96 ± 0.03% in Fig. 4(b); for the balanced 
functions, the contrast ratio 77 equals the interference visibility; in Fig. 4(c), 775 = 95.76 ± 0.07%; and in Fig. 4(c), 
i]b = 96.13 ± 0.07%. From Fig. 4(a) and 4(b), we can see that the photon counts of D2 fall with increasing voltage. 
This phenomenon is mainly caused by the coupling of multi mode fibers used in the detection part. We modulate 
the phase by changing the angle of the first mirror (changing the voltage of PZT). Although the change of the angle 
is very tiny, it will also affect the coupling efficiency, becoming worse when photons pass though the setup. Our 
experimental errors are mainly caused by the imperfections of PBS and HWP, the interference visibility, and the 
effect of DP [22I [231] - However, these errors can be reduced with improvement in the experimental technique. 

In conclusion, we have experimentally realized Deutsch's algorithm using linear optical components. We can 
determine a property of a function in one evaluation in the quantum case instead of two in the classical case. When 
phase ip = (2N + 1)tt, we need only a single photon as the input to judge the function f(x): a constant function 
when the photon is in V polarization and a balanced function when the photon is in H polarization. For the other 
input states, ip 7^ (2N + l)7r, we can only probably discriminate the function. We implement the CNOT gate using a 
Sagnac interferometer in the SPTQ logic. This experimental system is very stable and the experimental data are in 
excellent accordance with theoretical results. We believe these can be used to perform more complex entangled states 
or few-qubit quantum computation. 
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